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1 Introduction 



"There are more things in heaven and earth, Horatio, 
Than are dreamt of in your philosophy." 

"Hamlet", W. Shakespear 

In our first Chapter we are going to present some general notions and ideas of modern 
theory of integrability trying to outline its computability aspects. The reason of this approach 
is that though this theory was wide and deeply developed in the last few decades, its results 
are almost unusable for non-specialists in the area due to its complexity as well as due to some 
specific jargon unknown to mathematicians working in other areas. On the other hand, a lot 
of known results are completely algorithmical and can be used as a base for developing some 
symbolic programm package dealing with the problems of integrability. Creation of such a 
package will be of a great help not only at the stage of formulating of some new hypothesis but 
also as a tool to get new systematization and classification results, for instance, to get complete 
lists of integrable equations with given properties as it was done already for PDEs with known 
symmetries 

Some results presented here are quite simple and can be obtained by any student acquainted 
with the basics of calculus (in these cases direct derivation is given) while some ideas and results 
demand deep knowing of a great deal of modern mathematics (in this cases only formulations 
and references are given). Our main idea is not to present here the simplest subjects of inte- 
grability theory but to give its general description in simplest possible form in order to give 
a reader a feeling what has been done already and what could/should be done further in this 
area. Most of the subjects mentioned here will be discussed in details in the next Chapters. 

We will use the word "integrability" as a generalization of the notion "exactly solvable" 
for differential equations. Possible definition of differential operator will be discussed as well 
as some definitions of integrability itself. Numerous examples presented here are to show in 
particular that it is reasonable not only to use different notions of integrability for different 
differential equations but sometimes it proves to be very useful to regard one equation using 
various definitions of integrability, depending on what properties of the equation are under 
the study. Two classical approaches to classification of integrable equations - conservation 
laws and Lie symmetries - are also briefly presented and a few examples are given in order 
to demonstrate deep difference between these two notions, specially in case of PDEs. Two 
interesting semi-integrable systems are introduced showing one more aspect of integrability 
theory - some equations though not integrable in any strict sense, can be treated as "almost" 
integrable due to their intrinsic properties. 

2 Notion of differential operator 

There are many ways to define linear differential first-order operator D, starting with Leibnitz 
formula for product differentiating 
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This definition leads to 



/ 9 9-'" 

D ■ a = a + aD, D a = aD + 2a D + a , 



Dn ' a = Yl ( I ) ^W^"^ 



where 

n \ n(n — l)...(n — k + 1) . , / n 



with 1=1 



k J 1 • 2 • • • k V 

are binomial coefficients. Thus, one linear differential first-order operator and its powers are de- 
fined. Trying to define a composition of two linear operators, we get already quite cumbersome 
formula 

D v D 2 a = D x {D 2 (a) + aD 2 ) = D x {a)D 2 + D 1 D 2 {a) + D 2 (a)D 1 + aD x D 2 , (2) 

which can be regarded as a definition of a factorizable linear differential second-order op- 
erator. Notice that though each of D% and D 2 satisfies Leibnitz rule, their composition D\D 2 
does not! An important notion of commutator of two operators [D\D 2 — D 2 D\] plays a 
role of special multiplication 1 . In case of non-factorizable second-order operator an attempt 
to generalize Leibnitz rule leads to a very complicated Bourbaki-like constructions which we 
are not going to present here. All these problems appear due to coordinateness of this approach. 

On the other hand, a linear differential operator being written in coordinate form as 

j 

leads to 

L=^f a d a , d a = d ai ■■■d an , \a\ = orx + ... + a n 

\a\<m 

which is definition of LPDO of order m with n independent variables. Composition of two 
operators L and M is defined as 

L o M = f«d a 9p9 p = M 7 

and coefficients /i 7 are to be found from formula (JT]). 

Now, notion of linear differential equation (LDE) can be introduced in terms of the kernel 
of differential operator, i.e. 

Ker(L) := {tp\ ^ f a ff*<p = 0}. 

\a\<m 

Let us regard as illustrating example second-order LODO with one independent variable x: 

d 



L = f + f 1 d + f 2 d 2 , d:~- 

ax 



1 see Ex.1 
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Notice that if two functions ifii and ip 2 belong to its kernel Ker(L), then 

civ?i + c 2 (fi2 £ Ker(L), 
i.e. Ker(L) is a linear vector space over constants' field (normally, it is lorC). 



Main theorem about ODEs states the existence and uniqueness of the Cauchy problem for 
any ODE, i.e. one-to-one correspondence between elements of the kernel and initial data. In 
our case Cauchy data 



X=XQ 



dlp\ 



X=XQ 



form a two-dimensional vector space and, correspondingly, dimension of kernel is equal 2, 
dim(Ker(L)) = 2. Any two functions <pi,(p 2 £ Ker(L) form its basis if Wronskian < y?i,y?2 > 
is non-vanishing: 



¥1 ¥2 



7^0 



while an arbitrary function ip G Ker(L) has to satisfy following condition: 

< <pi,<p 2 , 4> >■- 



Pi ip 2 ip 
<p'i ¥2 

u 11 . // 
<Pl <P2 W 



0. 



Now we can construct immediately differential operator with a given kernel as 

L{$) =< <pi,<p 2 ,ij> > ■ 

For instance, if we are looking for an LODE with solutions sin a; and y/x, then corresponding 
LODE has form 

1 x 1 , 3 1 , 

ip (1 tanx) + tanxw w — -— ip — 0. 

2 2x Ax 2 

Coming back to LODO of order m, we re- write formula for the kernel in the form 



(3) 



< </?i, ...,(p m > 

which provides that high-order coefficient f m — 1. It is done just for our convenience and we 
will use this form further. 



Now, some constructive definition of linear differential operator was given and importance 
of its kernel was demonstrated. Corresponding differential equation was defined in terms of this 
kernel and construction of operator with a given kernel was described. All this is not possible 
for nonlinear operator because in this case manifold of solutions has much more complicated 
structure then just linear vector space - simply speaking, the reason of it is that in this case 
linear combination of solutions is not a solution anymore. Due to this reason some other notions 
are to be used to study properties of nonlinear operators - symmetries, conservation laws and, 
of course, as the very first step - change of variables transforming a nonlinear operator into a 
linear one. We will discuss all this in the next sections. 
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3 Notion of integrability 

• 3.1. Solution in elementary functions: 



y" + y = Q- 

General solution of this equation belongs to the class of trigonometrical functions, y = 
a sin(x + b), with arbitrary const a, b. In order to find this solution one has to notice that 
this equation is LODE with constant coefficients which possess fundamental system of 
solutions, all of the form e Xx where A is a root of characteristic polynomial. 

• 3.2. Solution modulo class of functions: 



3.2.1. 

y" = f(y)- 



In order to integrate this equation let us notice that y"y' = f(y)y' which leads to 

7^y' 2 = J f{y)dy + const = F(y) 



and finally 

dy 



dx 



which describes differential equation with separable variables. In case when F(y) 
is a polynomial of third or fourth degree, this is definition of elliptic integral 

and therefore the initial nonlinear ODE is integrable in elliptic functions. Particular 
case when polynomial F{y) has multiple roots might leads to a particular solution 
in elementary functions. Let us regard, for instance, equation y" = 2y 3 and put 
const = 0, then y' 2 = y 4 and y = 1/x, i.e. we have ONE solution in the class 
of rational functions. General solution is written out in terms of elliptic functions. 
Conclusion: equation is integrable in the class of elliptic functions and not inte- 
grable in the class of rational functions. 



Notice that as a first step in finding of solution, the order n of initial ODE was 
diminished to n — 1, in our case n = 2 . Of course, this is not possible for any 
arbitrary differential equation. This new ODE is called first integral or conserva- 
tion law due its physical meaning in applications, for instance, our example can be 
reformulated as Newton second law of mechanics and its first integral corresponds 
to energy conservation law. 



y" = y 2 + x. 

This equation defines first Painleve transcendent. About this equation it was 
proven that it has no solutions in classes of elementary or special functions. On 
the other hand, it is also proven that Painleve transcendent is a meromorphic 
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function with known special qualitative properties (fHj). 



This example demonstrates us the intrinsical difficulties when defining the notion 
of integrability. Scientific community has no general opinion about integrability of 
Painleve transcendent. Those standing on the classical positions think about it 
as about non-integrable equation. Those who are working on different applicative 
problems of theoretical physics involving the use of Painleve transcendent look 
at it as at some new special function . 

• 3.3. Solution modulo inexplicit function: 



u t = 2uu x . 

This equation describes so called shock waves and its solutions are expressed in terms of 
inexplicit function. Indeed, let us rewrite this equation in the new independent variables 
t = t, x = u and dependent one u = x, i.e. now x = 9(t, u) is a function on t, u. Then 



dt = dt, dx = 6 t dt + 9 u du, u t \d x =o = — 



0,, 



U X \dt=0 



and 



9 1 

= 2u— -6 t = 1u -9 = 2ut - cp(u) 

9n, 9n 



and finally x + 2tu = ip(u) where tp(u) is arbitrary function on u. Now, we have finite 
answer but no explicit form of dependence u = u(x,t). Has the general solution been 
found? The answer is that given some initial conditions, i.e. t = 0, we may define solution 
as u = ip~ l (x) where is tp~ l denotes inverse function for ip. 

3.4. Solution modulo change of variables (C-Integrability): 



0. 



This equation is called Thomas equation and it could be made linear with a change of 
variables. Indeed, let ip = log9 for some positively defined function 9: 



9j_ 
J 



1k 



ip x = (log0) x = — , ip y = (log9) y 
9 X y9 9y9 x 



<9 x . 



9y 

9' 

9 X 9y 



'xy \ q /y q 2 ' rxry q 2 

and substituting this into Thomas equation we get finally linear PDE 



9 xy + a9 x J t 



0. 



Suppose for simplicity that (3 = and make once more change of variables: 9 
then 

Ox = <Pxe klV } 9 xy = (j) xy e kiy + k 1( j) x e kiy and (j) xy + (h + a)(j) x = 0, 
and finally 



?kiy 







x\Vy 



+ (h + a)<p) = 0, 



which yields to 

4> y - k 2 <p = f(y) with k 2 = -(fa + a) 

and arbitrary function f(y). Now general solution can be obtained by the method of 
variation of a constant. As a first step let us solve homogeneous part of this equation, 
i.e. 4> y — k 2 4> = and 4>(x,y) = g(x)e k2V with arbitrary g(x). As a second step, suppose 
that g(x) is function on x, y, i.e. g(x,y), then initial equation takes form 

(g(x,y)e k *y) y -k 2 g(x,y)e k >y = f(y), 

g(x,y) y e k2V + k 2 g(x,y)e k2V - k 2 g(x,y)e k2y = f(y), 

g(x, y) y e k ™ = f(y), g(x, y) = J f^e^dy + h(x), 

and finally the general solution of Thomas equation with a = can be written out as 

<j>(x,y) = g(x,y)e k2 y = e k2 \ J f(y)e~ k2 My + h(x)) = f(y) + e^x) 

with two arbitrary functions f(y) and h(x). 
• 3.5. Solution modulo Fourier transform (F-Integrability): 

u t = 2uu x + eu xx 



where e is a constant. This equation is called Burgers equation and it differs from 
Thomas equation studied above where change of variables was local in a sense that solution 
in each point does not depend on the solution in some other points of definition domain, 
i.e. local in (x, y)-space. This equation can be transformed into 

u t = 2uu x + u xx 

by the change of variables x = ex, u = eu and t = e 2 t. We will use this form of Burgers 
equation skipping tildes in order to simplify the calculations below. To find solution of 
Burgers equation one has to use Fourier transform which obviously is nonlocal, i.e. here 
solution is local only in /c-space. In order to demonstrate it let us integrate it using 
notation 

h x=v - 

then after integration 

v t = v 2 x + v xx = e~ v (e v ) xx 
and change of variables w = e v , Burgers equation is reduced to the heat equation 

w t = w xx 

which is linear. Therefore, solutions of Burgers equation could be obtained from solutions 
of heat equation by the change of variables 

U! x 

U = V x = 

W 
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as 

, x f 00 . . /7\77 . . f exp(ikx - k 2 t)w(k)ikdk 

w(x, 0) = / exp(ikx)w(k)dk =>• i) = ^ V- , . . ' , . ,, 

v 7 v 7 w v y / exp(?£;x - k 2 t)w{k)dk 

where 0) is initial data and is a function called its Fourier transform and it 

can be computed as 



w 



1 f°° 

(k) = — / exp(— ikx)w(x, 0)dx. 

27T 



In fact, it is well-known that any linear PDE with constant coefficients on an infinite 
line can be solved using as standard basis {e tkx \k G R} because they are eigenfunctions 
of these operators. Thus, Thomas equation with a = where the general solution was 
found explicitly, is an exception while heat equation demonstrates the general situation. 



• 3.6. Solution modulo 1ST (S-Integrability): 

VLt = QztU x ~\~ Vb xxx . 

This equation is called Korteveg-de Vries (KdV) equation and it is nonlinear PDE 
with nonconstant coefficients. In this case, choosing set of functions {e tkx } as a basis is 
not helpful anymore: Fourier transform does not simplify the initial equation and only 
generates an infinite system of ODEs on Fourier coefficients. 



On the other hand, some new basis can be found which allows to reduce KdV with rapidly 
decreasing initial data, u — > 0, x — > ±oo, to the linear equation and to solve it. This new 
basis can be constructed using solutions of linear Schrodinger equation 



j>xx + k 2 jj = m\) 

where function u is called potential due to its origin in quantum mechanics. Solutions 
of linear Schrodinger equation are called Jost functions, ^(t^x^k), with asymptotic 
boundary conditions: 



^(t, x, k; u(x, t )) e - ±i{ - kx+kH) -> 1, x -> ±oo. 
Jost function (f>(x, k) = ip + (t, x, k; u(x, t))e~'^ kx+kit "> is defined by the integral equation 

k) = 1 + r 1 ~ eXp[ l k } X ~ X ' )] u(x') ( t ) (x > , k)dx' 
Jx 2/c 

with t playing role of a parameter. Second Jost function is defined analogously with 
integration over [— oo,x]. Notice that asymptotically for x — > ±oo linear Schrodinger 
equation 

i^xx + k 2 i[) = uip is reduced to ip xx + k 2 if) = 

as in case of Fourier basis {e lkx \k 6 1}. It means that asymptotically their solutions do 
coincide and, for instance, any solution of linear Schrodinger equation 

rj} ~ Cie i(kx+k : H) + C2e -i(kx+k : H) for x _^ 
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It turns out that solutions of KdV can be regarded as potentials of linear Schrodinger 
equation, i.e. following system of equations 

Ut = 6uu x + u xxx , 
Ipxx + k 2 ^ = wp 

is consistent and any solution of KdV can be written out as an expansion of Jost functions 
which in a sense are playing role of exponents e lhx in Fourier basis [TB*] . 

On the other hand, there exists a major difference between these two basis: Fourier basis 
is written out in explicit form via one function while Jost basis is written out in inexplicit 
form via two functions with different asymptotic properties on the different ends of a 
line. The crucial fact here is that two Jost functions are connected by simple algebraic 
equation: 

ip-{x, k, t) = a{k)ip + {x, —k, t) + b(k)e ik3t ip + (x, k, t) 

while it allows us to construct rational approximation of Jost functions for given a(k) 
and b(k) and, correspondingly, general solution of KdV. The problem of reconstruction 
of function u according to a(k),b(k) is called inverse scattering problem and this 
method, correspondingly, inverse scattering transform (1ST). 



• 3.6. Solution modulo Dressing method (D-Integrability): 

iu t = u xx ± \u\ 2 u. 



This equation is called nonlinear Schrodinger equation (NLS) and it is very important 
in many physical applications, for instance, in nonlinear optics. Dressing method is 
generalization of 1ST and in this case role of auxiliary linear equation (it was linear 
Schrodinger equation, second order ODE, in the previous case) plays a system of two 
linear first order ODEs [TK] : 

^ = + U ip® 
^ x 2) = -A^ (2) + W> (1) 

where v = ±-u. It turns out that system of equations 

{iut = u xx ± \u\ 2 u 
if,® = -Xtp® + vijjV 

is consistent and is equivalent to Riemann-Hilbert problem. Solutions of this last system 
are called matrix Jost functions and any solution of NLS can be written out as an 
expansion of matrix Jost functions which in a sense are playing role of exponents e lkx in 
Fourier basis ^2] • Detailed presentation of Dressing method will be given in some of our 
further Chapters. 
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4 Approach to classification 

Our list of definitions is neither full nor exhaustive, moreover one equation can be regarded 
as integrable due to a few different definitions of integrability. For instance, equation for 
shock waves from § 3.3, u t = 2uu x , is a particular form of Burgers equation from § 3.5, 
u t = 2uu x + eu xx with e — 0, and it can be linearized as above, i.e. it is not only integrable in 
terms of inexplicit function but also C-integrable and F-integrable, with general solution 

. . J exp(ikx — k 2 t)u(ek)ikdk 
J exp(ikx — k 2 t)u(ek)dk 

What form of integrability is chosen for some specific equation depends on what properties 
of it we are interested in. For instance, the answer in the form of inexplicit function shows 
immediately dependence of solution form on initial conditions - graphically presentation of 
inverse function u = ip~ 1 {x) can be obtained as mirrored image of x = <p(u). To get the same 
information from the formula above is a very nontrivial task. On the other hand, the general 
formula is the only known tool to study solutions with singularities. This shows that definitions 
of integrability do not suit to serve as a basis for classification of integrable systems and some 
more intrinsic ways should be used to classify and solve them. Below we present briefly two 
possible classification bases - conservation laws and Lie symmetries. 



4.1 Conservation laws 

Some strict and reasonable definition of a conservation law (which is also called first integral 
for ODEs) is not easy to give, even in case of ODEs. As most general definitions one might 
regard 

j f F(y) = for ODE j f y = f(y), y=( yi ,...,y n ), /(y) = (/i,...,/„) (4) 

and 
d f 

— G(u,u x ,u y ,u xx ,u xy ,u yy ...)dxdy... = for PDE d t u = g(u,u x ,Uy,u xx ,u X y,u yy , ...). 

(5) 

Obviously, without putting some restrictions on function F or G these definitions are too general 
and do not even point out some specific class of differential equations. For instance, let us take 
any second order ODE, due to well-known theorem on ODEs solutions we can write its general 
solution in a form 

F(t,y,a,b) = 

where a, b are two independent parameters (defined by initial conditions). Theorem on inexplicit 
function gives immediately 

a — Fi(t,y,b) Wb and b = F 2 (t,y,a) Va, 

i.e. any second order ODE has 2 independent conservation laws and obviously, by the same 
way n independent conservation laws can be constructed for ODE of order n. For instance, in 
the simplest case of second order ODE with constant coefficients, general solution and its first 
derivative have form 



y = cie 



+ c 2 e X2X 



CiAie AlX + c 2 \ 2 e X2X 
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and multiplying y by Ai and A2 we get equations on c\ and c 2 correspondingly: 



'\2y-y = c x (A 2 - A x )e Aia 
Aiy - 2/ = c 2 (Ai - A 2 )e A2a 



A 2 x + c 2 = - Xxy) 

\ lX + c"i = log(y - \ 2 y) 



and two conservation laws are written out explicitly: 

c 2 = log(y - Aiy) - A 2 x, c 1 = log(y - \ 2 y) - \ix. 

To find these conservation laws without knowing of solution is more complicated task then to 
solve equation itself and therefore they give no additional information about equation. This 
is the reason why often only polynomial or rational conservation laws are regarded - they are 
easier to find and mostly they describe qualitative properties of the equation which are very 
important for applications (conservation of energy, momentum, etc.) 

On the other hand, conservation laws, when known, are used for construction of ODEs 
solutions. Indeed, let us rewrite Eq.(jU) as 



d 
dt 



F(y) = (fidi + f 2 d 2 + ... + f n d n )F = C{F) = 0, 



i.e. as an equation in partial derivatives C(F) = 0. Such an equation has (n— 1) independent 
particular solutions ((p\, (p 2 , ip n -i) if its Jacobian matrix has maximal rank 



rank 



d{<Pi,<fi 2 , ...,<p n -i) 
d(yi,V2,...,y n -i) 



n 



(6) 



with notation 







/ &£± 




dtpi 






' 9yi 




dyn 


9(yi,y 2 ,- 


;Vn-l) 






dtfin 






\ dyi 


dyi 


dy n 



Now we can write out the general solution as F(ipi, <p 2 , 
and initial ODE can be reduced to 



tp n -i) with arbitrary function F 



dz dz 

Tt =f{zh W) 

and solved explicitly in quadratures (see § 3.2). 



dt 



4.2 Symmetry properties 

In order to give definition of canonical form for n-order ODE 

y^ = F(x,y,y',...,y^) 

let us first introduce vector 



(7) 



x 
V 

y 

Jn-l) 
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which is called vector of dynamical variables with all its coordinates regarded as indepen- 
dent, and its first derivative 



dt 



with respect to some new independent variable t such that dt = dx, then the equation 



" 1 " 




"1" 


y 




y 


y 




II 

y 


y( n ) 




_F_ 



dy 
dt 



1 

y 

ii 

y 



is called canonical form of an ODE. This canonical form is also called dynamical system, 

important fact is that dimension of dynamical system is (n + 1) for n-order ODE. 



Definition 4.1. Dynamical system 

f = 9{y) (9) 

dr 

is called a symmetry of another dynamical system 

dy 

dt = m 

d .dy. d .dy. d . d . ... . , 

holds. Symmetry g of dynamical system / is called trivial if g = const •/. 

Obviously, Eq. fjlUj) gives necessary condition of compatibility of this two dynamical sys- 
tems. It can be proven that this condition is also sufficient. Therefore, construction of each 
Lie symmetry with group parameter r is equivalent to a construction of some ODE which have 
t as independent variable and is consistent with a given ODE. System of ODEs obtained this 
way, when being written in canonical form is called dynamical system connected to the 
element of Lie symmetry algebra. 

We will discuss it in all details later (also see ^1]), now just pointing out the fact that Lie 
algebra can be generated not only by Lie transformation group (normally used for finding of so- 
lutions) but also by the set of dynamical systems (jHJ) (normally used for classification purposes). 

Two following elementary theorems show interesting interconnections between conservation 
laws and symmetries in case of ODEs. 
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Theorem 4.1 Let dynamical system 



dr 

is a symmetry of another dynamical system 

dy 
dt 



and 



d_ 

dt 



F{y) = (fid! + f 2 d 2 + ... + f n d n )F = C{F) = 



is a conservation law. Then Fg is symmetry as well. 
► Indeed, let us introduce 

M = (gid x + g 2 d 2 + ... + g n d n ), 
then consistency condition Eq.flTU]) is written out as 

Co M =M oC 

and after substituting FM. instead of M. we get on the left hand of (fTTj) 
C(FM) = C(F)M + FCM = FCM = FMC, 
and on the right hand of 

{FM)C = FMC. 



Corollary 4.2 Let n-th order ODE of the form has a symmetry 



dt 



go 

92 
9n 



then g Q is a conservation law and consequently without loss of generality we may put go — 1, if 

Theorem 4.3 An ODE of arbitrary order n having of (n — 1) independent conservation laws 
(i.e. complete set), has no nontrivial symmetries consistent with conservation laws. 

► Indeed, using all conservation laws we can reduce original n-order ODE into first order 
ODE 

da ... 
dt = /(a) 



and look for symmetries in the form 



da 
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Then 



and finally 



J ln(—) = const, 



f(a) 9(a) 9 
i.e. functions f(a) and (7(a) are proportional. I 

Simply speaking, Eq.(J5j) defines a one-parameter transformation group with parameter r 
which conserves the form of original equation. The very important achievement of Lie was 
his first theorem giving constructive procedure for obtaining such a group. It allowed him to 
classify integrable differential equations and to solve them. The simplest example of such a 
classification for second order ODEs with two symmetries is following: each of them can be 
transformed into one of the four types 

(I) y" = h(y') (II) y" = h(x), (III) y" = (IV) y" = h(x)y , 



where h denotes an arbitrary smooth function while explicit form of corresponding trans- 
formations was also written out by Lie. 

Symmetry approach can also be used for PDEs. For instance, first-order PDE for shock 
waves Ut = 2uu x (§ 3.3) has following Lie symmetry algebra u T = <p(u)u x where any smooth 
function (p = (p(u) defines one-parameter Lie symmetry group with corresponding choice of 
parameter r. Indeed, direct check gives immediately 

(u t ) T = (2(pu x )u x + 2u((pu xx + (p'ul) 

and 

(u T ) t = ip u x 2uu x + tp{2u 2 x + 2uu xx ) 

while relation {ut) T = (u T ) t yields to the final answer (here notation cp' = <p u was used). In 
order to construct dynamic system for this equation, let us introduce dynamical independent 
variables as 

u, u x = u x , u 2 = u xx , ... 

and dynamical system as 





u 




2uu\ 


d 


Ml 




2uu 2 + 2u\ 


dt 


u 2 




2uus + 6uiu 2 



This system can be transformed into finite-dimensional system using characteristics method 

|E|. ' 



For PDE of order n > 1 analogous dynamic system turns out to be always infinite and 
only particular solutions are to be constructed but no general solutions. Infinite-dimensional 
dynamical systems of this sort are not an easy treat and also choice of dynamical variables 
presents sometimes a special problem to be solved, therefore even in such an exhaustive textbook 
as Olver's ^T] these systems are not even discussed. On the other hand, practically all known 
results on classification of integrable nonlinear PDEs of two variables have been obtained using 



14 



this approach (in this context the notion of F-Integrability is used as it was done to integrate 
Burgers equation, § 3.5.) For instance, in ^B] for a PDE of the form 



u 



f(x,y) 



with arbitrary smooth function / on the right hand it was proven that this PDE is integrable 
and has symmetries iff right part has one of the following forms: e u , sinw or C\e u + C2e~ 2u with 
arbitrary constants C\,C2- Another interesting result was presented in jlO, where all PDEs of 
the form 

have been classified. Namely, a PDE of this form is integrable and has symmetries iff if it can 
be linearized by some special class of transformation. General form of transformation is written 
out explicitly. 



4.3 Examples 

Few examples presented here demonstrate different constellations of symmetries (SYM), con- 
servation laws (CL) and solutions (SOL) for a given equation(s). 

4.3.1: SYM +, CL +, SOL +. 

Let us regard a very simple equation 

y =i, 

then its dynamical system can be written out as 



dt 



"1" 




"i" 


y 




y 


F 




l 



with dt = dx 



and its general solution is y = ~x + c\x + c 2 with two constants of integration. In order 
to construct conservation laws, we need to resolve formula for solution with respect to the 
constants c\, c<i- 

ci = y -x, c 2 = y + -x 2 - xy' . 
Now, we look for solutions F(y) of the equation 

d 



dt 



(d x + yd y + d y/ )F = C(F) = 



with £ = d x + y'dy + d >. Direct check shows that F = y — x and F = y + \x 2 — xy' are 
functionally independent solutions of this equation. Moreover, general solution is an arbitrary 
function of two variables ^ 

F = F(y — x, y + -x 2 — xy'), 



for example, 



F = (y - x) 2 - 2(y + ^x 2 - xy') = y - 2y. 
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On the other hand, if there are no restriction on the function F, the conservation laws may 
take some quite complicated form, for instance, 

F = Arcsin(y — x)/(y+ -x 2 — xy') om . 

Now, that dynamical system, conservation laws and general solution of the original equation 
have been constructed, let us look for its symmetry: 




= g(y), g(y) = (91,92, gs)- 



Demand of compatibility 

= T^ 9 ^ is ec l uivalent to £(9i) = £(93) = °, £(92) = 93, 

and it can be proven that any linear combination of two vectors (1, 0, 0) and (0, x, 1) with (some) 
scalar coefficients provides solution of compatibility problem. Thus, Lie symmetry group cor- 
responding to the vector (1,0,0) is shift in x while the second vector (0,x, 1) corresponds to 
summing y with particular solutions of homogeneous equation. 

4.3.2: SYM -, CL +, SOL +. 

Let us regard system of ODEs 

= (n 2 - n 3 )a 2 a 3 

= (n 3 - ni)aia 3 (12) 
= (m - n 2 )a 1 a 2 

with variables dj and constants n*. This system is well-known in physical applications - it 
describes three-wave interactions of atmospheric planetary waves, dynamics of elastic pendulum 
or swinging string, etc. Conservation laws for this system are: 

• Energy conservation law is obtained by multiplying the i-ih equation by a i? i = 1,2,3 
and adding all three of them: 

ilia 2 + n2d2 + n^a\ = const . 

• Enstrophy conservation law is obtained by multiplying the i-th equation by rijOj, 2 = 1,2,3 
and adding all three of them: 

n\a\ + n\a\ + n\a\ = const . 

Using these two conservation laws one can easily obtain expressions for a 2 and a 3 in terms 
of a\. Substitution of these expressions into the first equation of Svs. (fT2^) gives us differential 
equation on ai: 

(y'f = f(y) whh y = ai 

whose explicit solution is one of Jacobian elliptic functions while a 2 and a 3 are two other Ja- 
cobian elliptic functions jH]. In fact, Svs.(fl2"|) is often regarded as one of possible definitions of 
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Jacobian elliptic functions. 



Notice that from Theorem 4.3 one can conclude immediately that all symmetries consis- 
tent with conservation laws, are trivial. 



4.3.3: SYM +, CL ±, SOL +. 

Let us regard heat equation 

d t u = d xx u, (13) 

which generates solutions of Burgers equation (see § 3.5). Direct check shows that Eq. (fT3*j) 

is invariant due to transformations x = xt and t = tr 2 with any constant r, i.e. dilation 

transformations constitute Lie symmetry group for Eq. ljlHj) . Moreover, Eq. ljlHj) is integrable 
and its solution 



u 



/exp(ikx — k 2 t)u(k)dk with u(k) = — / exp(— ikx)u(x, 0)dx 
27T 



is obtained by Fourier transformation. 

This example demonstrates also some very peculiar property - heat equation (as well as 
Burgers equation) has only one conservation law: 

d f°° 

— / udx = 0. 

J— oo 

Nonexistence of any other conservation laws is proven, for instance, in [1 



5 Semi-integrability 

5.1 Elements of integrability 

Dispersive evolution PDEs on compacts is our subject in this subsection. In contrast to stan- 
dard mathematical classification of LPDO into hyperbolic, parabolic and elliptic operators 
there exists some other classification - into dispersive and non- dispersive operators - which is 
successfully used in theoretical physics and is not complementary to mathematical one (for 
details see |H]). Let regard LPDE with constant coefficients in a form 

where t is time variable and x is space variable, and suppose that a linear wave 

ijj(x) = A exp i(kx — ut) 

with constant amplitude A, wave number k and frequency u is its solution. After substituting a 
linear wave into initial LPDE we get P(—iu>,ik) = 0, which means that k and uj are connected 
in some way: there exist some function / such, that /(o>, k) = 0. 

This connection is called dispersion relation and solution of the dispersion relation is 
called dispersion function, u — u(k). If condition 
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d 2 u 



7^0 



dk 2 



holds, then initial LPDE is called evolution dispersive equation and it obviously is com- 
pletely defined by dispersive function. All these definitions can be easily reformulated for a 
case of more space variables, namely Xi,X2, in this case linear wave takes form 



dispersion function can be computed from P{— iui, iki, ik n ) = and the condition of 
non-zero second derivative of the dispersion function takes a matrix form: 



Notice now that solutions of linear evolution dispersive PDE are known by definition and 
the reasonable question here is: what can be found about solutions of nonlinear PDE 



with dispersive linear part C(ip) and some nonlinearity J\f(ijj)7 

Nonlinear PDEs of this form play major role in the theory of wave turbulence and in general 
there is no final answer to this question. Case of weak turbulence, i.e. when nonlinearity 
jV(V') is regarded small in a sense that wave amplitudes A are small enough (smallness of an 
amplitude can be strictly defined), is investigated in much more details. Two qualitatively 
different cases have to be regarded: 

— * 

1. coordinates of wave vector are real numbers, {k = (k%, ....k n )\ki e R} (corresponds to 
infinite space domain); 

2. coordinates of wave vector are integer numbers, {k = (ki, ....k n )\ki G Z} (corresponds 
to compact space domain). 

In the first case method of wave kinetic equation has been developed in 60-th (see, for 
instance, [3]) and applied for many different types evolution PDEs. Kinetic equation is approx- 
imately equivalent to initial nonlinear PDE but has more simple form allowing direct numerical 
computations of each wave amplitudes in a given domain of wave spectrum. Wave kinetic 
equation is an averaged equation imposed on a certain set of correlation functions and it is 
in fact one limiting case of the quantum Bose-Einstein equation while the Boltzman kinetic 
equation is its other limit. Some statistical assumptions have been used in order to obtain 
kinetic equations and limit of its applicability then is a very complicated problem which should 
be solved separately for each specific equation [Tlj . 




with wave vector k = (k%, ....k n ) and space-like variable x 



(xi,. ..,£„). Then 




dkidkj 
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In the second case, exact solutions in terms of elliptic functions have been found [Zj . More 
precisely, it is proven that solving of initial nonlinear PDE can be reduced to solving a few 
small systems of ODEs of the form 





= £iA 2 Az 


A 2 




M 


= ^A 1 A 2 



(14) 



in case of quadratic nonlinearity, 



Mi 

A 2 
A 3 
A 4 



d x A 2 AsA A , 
<L 2 A X A Z A^ 
UA X A 2 A A , 
UA X A 2 A^ 



in case of cubic nonlinearity and so on. Notice, that in contrast to a linear wave with a constant 
amplitude A 7^ A(t, x), waves in nonlinear PDE have amplitudes A, depending on time. It 
means that solutions of initial nonlinear PDE have characteristic wave form as in linear case 
but wave amplitudes are Jacobian elliptic functions on time, cn(T),dn(T) and sn(T). Notice 
that Sys.(HU) has been studied in § 4.3.2 and its conservation laws were found. Exact solutions 
of Sys.(HU) are 

' A x = bicn(T/t - A), 
A 2 = b 2 dn(T/t - A), 
k ^3 = b 3 sn(T/t - A), 

and constants bi,t ,X are written out explicitly as functions of initial values of waves' ampli- 
tudes (see jH] for details). 



These systems of ODEs providing exact solutions of initial nonlinear evolution PDE are 
called elements of integrability. Some constructive procedure, Clipping method, has 
been developed (HI allowing to find all elements of integrability for a given evolution PDE. 



5.2 Levels of integrability 

Let us formulate classical three-body problem whose integrability attracted attention of many 
investigators beginning with Lagrange. Computing the mutual gravitational interaction of three 
masses is surprisingly difficult to solve and only two integrable cases were found. For simplicity 
we regard three-body problem with all masses equal, then equations of motion take form 



= ^12/12 + 213/13 
^Jr = ^21/12 + ^23/23 (15) 

= 231/13 + 232/23 

where Zj is a complex number, Zj = Xj + iyj, describing coordinates of j-th mass on a plane, 
fjk is a given function depending on the distance between j-th and fc-th masses (physically it 
is attraction force) while following notations are used: Zjk = Zj — z^ and fjk = f(\zjk\ 2 )- 



This system admits following conservation laws: 
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Velocity of center of masses is constant 

Summing up all three equations, we get 

d 2 

— Oi + z 2 + z 3 ) = 0. 

This equality allows us to choose the origin of coordinate system in such a way that 

z\ + z 2 + z 3 = 

which simplifies all further calculations significantly. That is the reason why till the end 
of this section this coordinate system is used. Physically it means that coordinate system 
is connected with masses center. 

Conservation of energy 

Multiplying j-th equation by z-, summing up all three equations and adding complex 
conjugate, we obtain on the left 

3 ,3 

i=i j=i 
i.e. left hand describes derivative of kinetic energy. 

On the right we have derivative of potential energy U: 

jU = j t (F(\z 12 \ 2 ) + F(\z 13 \ 2 ) + F(\z 23 \ 2 )) with notation F' = f 
and finally energy conservation law takes form 

3 

z'jz'j = F{\z 12 \ 2 ) + F{\z 13 \ 2 ) + F{\z 23 \ 2 ) + const 

3=1 



Conservation of angular momentum 

By differentiating of angular momentum 

3 



with respect to t, we get 



Im ( \ Zj\ 2 + f\ 2 Z\ 2 Z\ 2 + f\ 3 Z\ 3 Z\ 3 + ^23^23^23 

Im (/i2kl2| 2 + /l3 kl3 1 2 + /23k23| 2 ) = 0, 



while force / is some real-valued function. 

In general case there are no other conservation laws and the problem is not integrable. 
On the other hand, one may look for some periodical solutions of Svs.(fT5j) and try to deduce 
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the necessary conditions of periodicity. Importance of the existence of periodical solutions was 
pointed out already by Poincare and is sometimes even regarded as a definition of integrability 
- just as opposite case for a chaos. 

Theorem 5.1. If fij > 0, Vi, j (so-called repulsive case), then Svs. 4i<5j) has no periodical 
solutions. 

► Indeed, in case of periodical solution, magnitude of inertia momentum 

CT I 1 2 I I 1 2 I I 1 2 

'■— \Zl2\ + |^13 1 + | #23 1 

should have minimums and maximums as sum of distances between three masses. On the other 
hand, 

1 d 2 

(\ 1 2 I I 1 2 I I |2\ I • |2 | I • |2 i I • |2 i f I |2 i j? I |2 , f I |2 ^ n 
2,0^ 12 '^ 13 ' ' Z23 ' ! ~ \ Zl2 \ ~"~ I ^13 1 + |^23 1 +Jl2pl2| +Jl3pl3| +J23p23| > U, 

which contradicts to the fact that function Z has to have different signs in the points of mini- 
mum and maximum. I 

One interesting case - Poincare case - though does not lead to integrable reduction of 
Sys.(JT3J), give quite enlightening results and allows to regard this case as "almost" integrable. 
In this case there exists one more conservation law - conservation of inertia momentum 

1-212 | 2 + |-2l3 | 2 + 1-223 | 2 = Const 

and it is possible due to a special choice of function fjk = l/\zjk\ 4 which allows us to reduce 
initial system to the ODE of the form 

B" =a(B') 3 + b(B') 2 + cB' +d. (16) 

in new polar coordinates. This equation describes a geometrical place of points on the plane, 
i.e. some plane curve B, providing solutions of initial system. The curve B is of figure-eight 
form and can not be described by any known algebraic curve. On the other hand, it can be 
approximated with desirable accuracy, for instance, by lemniscate 

4 i 2 2 i n 4 2 2 

x + ax y + py = x — y . 
Very comprehensive collection of results and graphics one can find in ||] 

There exists hypothesis that the only periodical solution of Eq. (jl6J) is this eight-form 
curve (not proven). Existence theorem for non-equidistant periodical solutions is proven for a 
wide class of functions / (in variational setting). 

The simplest possible case of periodical solution can be obtained if one of Zj is equal to 
zero 2 (obviously, the problem is reduced to two-body problem). Two more complicated classical 
integrable cases with periodical solutions for particular choice of z-s are known: 

2 see Ex.6 
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Lagrange case: \zi 2 \ = \zrs\ = l^l- 

It means that distances between three masses are equal as well as all corresponding 
attraction forces, and the masses are moving along a circle. Lagrange case is also called 
equidistance case. In this case Svs.()15p can be reduced to ODE y = f(y) and solved in 
quadratures. 

Calogero case: all Zj are real and fjk = l/\zjk\ 4 for j, k = 1, ...n 

It means that all masses are moving along a line (in fact, along a real axes) and this 
is generalization of Euler case of three-body problem which after appropriate change of 
variables Svs. (jl5j) takes form 

d 2 xi _ 1 _|_ 1 
dt 2 (xi—x 2 ) 3 (xi—xs) 3 

d 2 x 2 _ 1 _|_ 1 



dt 2 (xi— x 2 ) 3 (x2-x 3 ) 3 

d 2 x 3 _ 1 1 

dt 2 (si— X3) 3 {X2-X3) 



Euler 's system was generalized in jT], [2] to 



with U = — and j — l,...,n. 



dt 2 dxj ^— ' (xi — x 

All n independent conservation laws were found and it was proven that the system is 
integrable. 



6 Summary 

In our first Chapter we introduced a notion of differential operator and gave few different def- 
initions of integrable differential equation. It was shown that some of them can be equivalent 
for a given equation and it is reasonable to choose an appropriate one depending on what prop- 
erties of the equation are under the study. Some interesting and physically important examples 
of "almost" integrable systems were described. Very intrinsic question on interconnections of 
conservation laws and symmetries was also discussed. 

Simplest possible example of nonlinear equation - famous Riccati equation which is ODE of 
first order with quadratic nonlinearity. We will use this equation in our next Chapter, Chap- 
ter 2: Riccati equation, in order to demonstrate many properties of differential equations 
described above. Riccati equation will also be very useful for introduction of some new notions 
like singularities of solutions, integrability tests, etc. 



7 Exercises for Chapter 1 

1. Using (J2J) prove that D3 = -D2-D1 — D1D2 satisfies Leibnitz rule. 

2. Prove that operator 

C = J2 a k* k d k x 

can be transformed into an operator with constant coefficients by the following change of 
variables: x = e*. (Euler) 



22 



3. Transform equation u% = ip(u)u x into Vt = vv x by appropriate change of variables. 

4. Prove that LODE with constant coefficients has conservation laws of the form 

l^lle-^ = const . 

5. For potential energy U from § 5.2 prove that (Lagrange- Jacobi identity) 

dU dU dU dU dU dU n/J! . l2 , . l2 r . l2 , 

Xl-K \-yi^ hX 2^ hX 3 - hj/3-a— =2(/i2 Z12 + /13 «13 + 723 ^23 )■ 

dxi dyi ox 2 oy 2 ox 3 dy 3 

6. For particular case z\ — and our choice of coordinate system solve Svs. (jl5|) explicitly. 
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